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SUMMAEY 


The  errors  involved  in  using  the  perfeci>-gaB  law  pv  =  ET  and  the 
assumption  of  constant  heat  capacities  are  evaluated.  The  "basic  equations 
of  gas  flows  taking  into  account  these  phenomena  separately  and  at  the 
same  time  are  presented. 


INTRODUCTION 


The  conventional  method  of  obtaining  high  Mach  numbers  for 
aerodynamic  tests  is  to  accelerate  the  air  by  means  of  a  pressure 
difference  so  that  the  random  kinetic  energy  of  the  molecules  of  air  at 
rest  is  converted  into  kinetic  energy  in  the  test  section.  For  very 
high  Mach  numbers  this  may  occasion  high  stagnation  temperatures  and 
pressures  which  introduce  effects  due  to  the  vibrational  heat  capacity 
and  molecular  forces  and  size  such  that  the  perfect— gas  law  pv  =  ET 
and  the  assumption  of  constant  heat  capacities  may  be  no  longer 
sufficiently  accurate  to  evaluate  gas  flows. 

It  is  the  purpose  of  this  paper  to  present  formulas  which  are 
suitable  for  handling  such  problems  and  to  point  out  the  magnitude  of 
the  errors  that  may  be  involved  in  using  the  perfect— gas  law  and  the 
assumption  of  constant  heat  capacities. 

Tsien  (reference  l)  has  published  a  theoretical  discussion  of 
this  problem  in  which  certain  approximations  were  introduced  in  order 
to  obtain  solutions  that  were  in  a  very  neat  form  when  the  imperfect- 
gas  effects  were  moderate.  A  comparison  of  Tsien' s  results  with 
this  work  is  presented  to  show  the  magnitude  of  these  approximations. 

In  England  Goldstein  has  previously  investigated  this  problem 
at  moderate  temperatures  and  pressures  in  order  to  prove  the  small 
magnitude  of  imperfect-gas  and  vibrational— heat— capacity  effects  in 
most  supersonic  wind  tunnels.  This  report  indicates,  in  general,  the 
range  in  which  these  effects  are  small  but  does  not  present  formulas  for 
handling  problems  in  gas  dynamics  when  these  effects  are  large. 
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The  present  paper  la  arranged  in  the  following  three  parts: 
temperature  effects  on  perfect— gas  flows  due  to  variation  of  heat 
capacities;  imperfect— gas  effects  on  gases  without  variation  of 
heat  capacities;  and  gas  flows  in  which  both  effects  are  present. 
This  is  done  since  the  formulas  in  the  first  part  may  prove  use- 
ful to  those  dealing  with  the  flow  of  hot  exhaust  gases  and  since 
it  may  bring  out  more  clearly  the  differences  between  the  two  effects. 


SYMBOLS 

a  term  in  Van  der  Waals'  equation  correcting  for  the 
effect  of  molecular  forces 

b  term  in  Van  der  Waals'  equation  correcting  for  the 

effect  of  molecular  size 

c  speed  of  sound,  feet  per  second 

Cp  heat  capacity  at  constant  pressure 

Cv  heat  capacity  at  constant  volume 

E  .  energy,  foot-pounds 

M  Mach  number   (w/c) 

p  pressure,  pounds  per  square  foot 

B  gas  constant 

T  absolute  temperature,  degrees  Fahrenheit 

v  specific  volume,  cubic  feet  .per  slug 

w  velocity,  feet  per  second 

y  ratio  of  heat  capacities  (Cp/CT) 

p  density,  slugs  per  cubic  foot 

6  characteristic  temperature  of  molecular  vibration 

Subscripts: 

o  stagnation  conditions 

c  critical  conditions 
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Errors  Involved  in  Assuming  Constant  Specific  Heats  in  the 

Presence  of  High  Temperatures  in  a  Perfect  Gas 

For  a  perfect  gas  with  constant  heat  capacities  the  equation  for 
conservation  of  energy  of  a  steady  isentropic  process  may  "be  written  as 


CpT  +  ^2  =  c?To 


If  this  equation  is  combined  with  the  equation  for  the  isentropic  speed 
of  sound  <fi   =  7BT  the  resulting  equation  is 


M^  = 


Tr 


7  -   1  \T 


-  1 


(1) 


If  the  expansion  is  isentropic,  the  pressure  and  density  ratios 
corresponding  to  the  Mach  number  are 


-fc 


and 


12 
P 


Cp/R 


T  Npv/k 


(2) 


However,  if  the  temperature  of  the  gas  is  high  enough  the  h9at 
capacities  may  not  he  assumed  constant  "because  the  vibrational  degrees 
of  freedom  of  polyatomic  molecules  are  excited.  The  variation  of  the 
equilibrium  value  of  the  heat  capacity  at  constant  volume  of  a  perfect 
diatomic  gas  is  found  from  quantum  mechanical  considerations  to  be  of 
the  form 


E 


5  + 
2 


0\2 


fl/T 


>A  -  1)' 


(3) 


where  6     is  a  constant  depending  on  the  gas.  The  formula  may  be  used 
for  the  mixture  air  if  the  value  of  9      is  placed  equal  to  5526  when 
absolute  temperature  is  measured  in  degrees  Fahrenheit.   (See  the 
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appendix.)  The  value  of  the  heat  capacity  at  constant  pressure  for  a 
perfect  gas  is  then 


*   2   W  (ee/T  _  l}2  M 


Figure  1  is  a  plot  of  equation  (k)   and  shows  that  the  heat  capacity 
may  not  he  considered  constant  above  600°  F  absolute. 

When  the  heat  capacity  at  constant  pressure  varies  according  to 
equation  (k)   the  energy  equation  must  he  written 


flT  c 
B  /   IE  dT  +  1*2  =  0  (5) 

UTo 
Substituting  equation  (k)    into  equation  (5)  and  integrating  yields 


^T  +  eeW-     +w2  =  'mo+    e/f°  (6) 


The  Mach  numher  is  obtained  from  equation  (6)  by  dividing  through 
by  7FT  =  c^}   which  gives 


e-lfi?-l)  + 


±  I  ^  -.  1  ]  +  f|  (_1 i )  (7) 


where 


The  pressure  ratio  corresponding  to  this  Mach  number  is  obtained  from 
the  iuentropic  equation 
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loSX  =    /T£P  AT  ,9) 

Po       jTo    E     T  ^ 

"by  substituting  equation   (Ij-)    into  equation   (9)   and  integrating  to  give 

e     e0/T        e      e0/Tc 
X  =  fX\7/2  1  -  e0/To  eVT  e0/T-1_  T°  eQ^-lJ 

Po     vW      x  _  ee/T 

Similarly,  the  density  ratio  is  found  to  "be 

ee/T      e     ee/TQ 


(10) 


p         /  T\5/2  x  _  eeAo     \T  eS/T-1     T0  ee/T0_^ 


(11) 


The  differences  involved  in  the  use  of  equations  (l)  and  (2)  to  predict 
the  temperature,  density,  and  pressure  ratios  corresponding  to  a  given 
Mach  number  are  given  in  figures  2(a)  and  2(h),  in  terms  of  the  percentage 
differences  from  the  value  given  "by  equations  (7),  (10),  and  (ll)  for 
stagnation  temperatures  of  1000°  and  2000°  F  absolute.. 

It  is  seen  that  the  assumption  of  constant  heat  capacity  leads  to 
appreciable  differences  in  applying  the  isentropic  law  for  a  perfect  gas 
if  stagnation  temperatures  above  1000°  F  absolute  are  involved. 


Errors  Involved  in  the  Assumption  of  the  Perfect-Gas 

Law  pv  =  ET  for  a  Gas  with  Constant  Heat  Capacities 

In  order  to. evaluate  flows  in  which  imperfect— gas  effects  are  present, 
an  equation  of  state  that  takes  into  account  these  effects  must  be  chosen. 
For  the  purposes  of  this  paper  an  equation  which  takes  into  account  the 
effects  of  molecular  forces  and  size  should  be  sufficient. 

A  suitable  equation  is  that  of  Van  der  WaalB 


P  +  -\\(r  -  b)  =  ET  (12) 
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where  b  is  a  term  correcting  for  the  volume  occupied  ty  the  molecules 
and  a  is  a  term  correcting  for  the  effect  of  molecular  forces. 

Figure  3  is  a  graph  of  Van  der  Waals'  equation  in  which  the  quanti- 
ties Pj  v ,   and  T  have  been  made  nondimensional  by  dividing  by  the  values 
of  these  quantities  at  the  critical  point  pCj  vc,  and  TCf   thus  making 
the  graph  suitable  for  any  gas.   (See  reference  2.)  The  graph  may  be  used 
for  air  if  an  empirical  critical  point  (pc  =  37.2  atm,  Tc  =  238. 5°  F  abs., 
vc  =  0.6^38  slugs /ft 3)  is  assigned  to  that  mixture  of  oxygen  and  nitrogen. 
To  give  this  critical  point  the  values  of  a,  b,  and  E  for  air  when  the 
pressure  is  measured  in  pounds  per  square  foot,  the  specific  volume  in 
cubic  feet  per  slug  and  the  absolute  temperature  in  degrees  Fahrenheit 
are  a  =  8.78  x  10-5,  b  =  0.65^,  and  E  =  1716. 

The  proper  equation  for  an  isentropic  expansion  of  a  real  gas  is 
(see  reference  3) 


dE  =  Cv  dT  + 


\dT/v=Constant 


dv  =  0 


(13) 


which  for  Van  der  Waals1  equation  becomes 


dE  =  Cv  dT  +  p  dv  +  -%  dv  =  0 

•7d 


(Ik) 


Equation  (lU)  may  be  written  as 


dE  =  Cv  dT  +  d(pv)  -  v  dp  +  -%  dv  =  0 

v^ 


and  since  — v  dp  =  w  dw 


dE  =  Cv  dT  +  d(pv)  +  -%  dv  +  w  dw  =  0 


(15) 


Assuming  constant  heat  capacity  at  constant  volume  and  integrating 
equation  (15)  gives 


w2 


E  =  CVT  +  v  (v   -  -%A  +  £■  =  Constant  =  EQ 


(16) 
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This  ie  then  the  energy  equation  for  a  Tan  der  Waals  gas.  Dividing 
through  by  the  isentropic  speed  of  sound 

c2  =  &  =    (l  +  £-\       t2bT ^  (17) 

dp   ^   CY)   (v  _  t)2   v  ^ 

and  since 

p  _  _&_  -  RT   _  2a 
t^   v  —  b   v2 


then 


VTo  ♦  ^o  (^  -  %)  -  2CTT  -  Br  (j^  -  £) 


2C 
M2 


The  value  of  v  for  an  isentropic  expansion  to  he  placed  in 
equation  (l8)  can  he  formed  from  equation  (l4)  as  follows: 


(V  dT  +  p  dv  +  -i  dv  =  (^  dT  +  R^  .  dv  =  0 


v 


then 


ii!  =  _jr_  (19, 

T      v  -  h 


and  if  Cv  is  constant 


/T  \°v^ 
v=(v0-h)(4pj     +  b  (20) 


From  equations  (l8)  and  (20) f   knowing  the  stagnation  conditions  for 
an  expansion  from  T0  to  T^  the  Mach  number  may  he  calculated.  The 
pressure  ratio  is  then  found  to  he 
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p 

RT             a 
t-Tj       t2 

Po 

ETo 

(21) 


v„  -  to 


and  substituting  the  value  of  v  from  equation  (20)  we  obtain 


Po 


(22) 


Figure  h   shows  the  conventional  pressure  ratio  and  area 
ratio  pw/(pw)j4_;L  plotted  against  Mach  number  for  air  starting  from 
stagnation  conditions  of  520°  F  absolute  and  various  pressures  compared 
with  the  value  obtained  using  constant  ratio  of  heat  capacities  and  the 
perfect— gas  law. 

Also  shown  in  figure  k   are  the  values  of  pw/(pw)j^-j_  computed  by 
Tsien's  method.   It  is  seen  that  as  the  imperfect— gas  effects  become 
large  it  is  no  longer  possible  to  simplify  the  analysis  by  neglecting 

terms  containing  the  squares  of  —  and  — ^—,   although  Tsien's  results 

v  pv 

are  in  good  agreement  at  50  atmospheres  when  the  Van  der  Waals  effect 
is  moderate. 

It  is  interesting  to  note  that  the  speed  of  sound  in  a  Tan  der  Waals 
gas 


dp    I1   Cj  (v  _  t)2    v  (17) 


is  not  equal  to  7FT.  The  expression  for  the  ratio  of  specific  heats  in 
a  Van  der  Waals  gas  is 

r-i  +  ff  r2  +  %ab^  Ce3) 

W  \  pT2  _  a  +  £M 
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The  value  of  the  ratio  of  heat  capacities  7     in  this  case  is 


2  +  (  T 


7  = 


6/T 


pv^  +  a 


(e^-lf   P-2~a  + 


2ab 

T 


l  +  M 


3e/T 


T^  (e*A-l)' 


(29) 


"but  the  speed  of  sound  is  found  from  equation  (17)  by  substituting  the 
value  of  Cv/R  from  equation  (3)  to  he 


C2  = 


1  + 


2  +  (  T 


.e/T 


(e*/T  -  l)' 


v^T 


(v  -  h)2 


2a 
v 


Figure  6  shows  the  conventional  pressure  ratio  and  area  ratio  plotted 
against  Maeh  number  for  air  starting  from  stagnation  conditions  of 
2000°  F  ahsolute  and  various  pressures  compared  with  the  value  obtained 
using  constant  ratio  of  heat  capacities  and  the  perfect— gas  law. 


DISCUSSION 


The  foregoing  analyses  show  that  the  effects  of  variation  of  heat 
capacities  with  temperature  do  not  "become  important  in  isentropic  expansions 
of  air  until  stagnation  temperatures  of  the  order  of  1000°  F  ahsolute  are 
encountered.  Above  1000°  F  ahsolute,  however,  for  accurate  analysis  this 
variation  must  he  taken  into  account.  In  general,  it  may  "be  stated  that 


for  diatomic  gases  these  effects  are  important  when  (— 
becomes  appreciable  compared  to  the  number  2.5. 


e/T 


'e/T 


-  1 


The  effects  of  "Van  der  Waals1  forces  become  important  when  either 
the  temperature  is  extremely  low  for  near  atmospheric  pressures  or  the 
pressure  very  high  for  moderate  temperatures.  These  forces  must  be 
taken  into  account  when  the  value  of  a/v2  becomes  appreciable  compared 
to  the  pressure  p,  or  b  becomes  appreciable  compared  to  v.  For  air 
these  effects  are  unimportant  until  stagnation  pressures  of  the  order 
of  50  atmospheres  at  stagnation  temperature  of  520°  F  absolute  are 
encountered. 
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Tsien's  method  agrees  well  with  the  results  of  this  investigation 
up  to  50  atmospheres  in  this  case,  hut  it  appears  that  it  is  not 

possible  to  neglect  the  squared  terms  of  —  and  — j|—  when  the  effects 


of  Van  der  Waals'  forces  become  appreciable. 


CONCLUSIONS 


v^p 


In  many  cases  found  in  very  high  Mach  number  wind  tunnels  and  in 
flows  of  high  stagnation  temperature  or  pressure,  imperfect— gas  effects 
and  the  effects  of  variation  of  heat  capacities  may  be  present. 

For  diatomic  gases  the  effect  of  variation  of  heat  capacities  becomes 

/fl\2    p0/t  / 

important  when  (  — )  s ^    becomes  appreciable  compared  to  5/2. 


For  air  these  effects  become  appreciable  when  stagnation  conditions  of 
1000°  F  absolute  or  larger  are  encountered. 

Imperfect— gas  effects  become  important  in  gas  dynamics  when  a/v2 
becomes  appreciable  compared  to  the  pressure  p  or  b  becomes  appreciable 
compared  to  v.  When  air  is  expanded  from  a  stagnation  temperature 
of  520°  F  absolute  these  effects  become  important  if  the  stagnation  pressures 
are  of  the  order  of  50  atmospheres  or  greater. 

Formulas  are  presented  for  handling  isentropic  expansions  taking  into 
account  these  phenomena  both  separately  and  at  the  same  time.  Tsien's 
method  is  found  to  be  applicable  for  small  departures  from  a  perfect  gas 
but  is  not  accurate  when  the  effects  of  Tan  der  Waals'  forces  become 
appreciable. 
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APPENDIX 
DERIVATION  OF  THE  VIBRATIONAL  HEAT  CAPACITY  OF  A  DIATOMIC  GAS 

To  arrive  at  the  vibrational  heat  capacity  of  a  diatomic  gas,  the 
individual  molecules  are  treated  as  linear  harmonic  oscillators  of  a 
fundamental  frequency  and  Shrodinger's  equation  is  solved  for  the 
allowable  energy  states  of  such  an  oscillator.  These  allowable  states 
are  then  substituted  into  the  equation  for  the  canonical  energy  distribu- 
tion and  the  average  energy  per  particle  as  a  function  of  the  absolute 
temperature  is  found.  This  may  be  differentiated  to  obtain  the 
contribution  of  the  vibrational  degrees  of  freedom  of  the  molecule  to 
the  heat  capacity  of  the  gas  at  any  temperature. 

The  average  vibrational  energy  per  particle  found  in  this  way  is 
(see  references  k   and  5) 


!  =  £¥-  +     hv 

2    ehv/kT.-L 

where 

h  Planck's  constant 

v  characteristic  frequency  of  molecular  vibration 

T  absolute  temperature 

Differentiating  to  obtain  the  contribution  to  the  heat  capacity  of  this 
energy  yields 


Cylb   1  M   /hvf    ehV^ 
E     k  ST  {yiJ     ^hv/kT  _  *2 

hv 
For  a  particular  gas  -s—  =  6     is  a  constant  and  may  be  determined  from 

spectroscopic  data.  The  heat  capacity  at  constant  pressure  is  then 


£e  -  1  +  Cyj-b  _  2  +  /£\2    ee/T 

E    2    E     2   \TJ     (    e/T   ,\< 

)   I     -1) 
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The  value  of  0  for  oxygen  iB  U010.U  and  for  nitrogen  is  60kh.h   for 
absolute  temperatures  measured  in  degrees  Fahrenheit.  The  value  5526 
may  "be  used  for  air. 
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Figure  2.-  Percent  error  involved  in  the  use  of  constant -heat-capacity 
formulas  to  obtain  T/T  ,   p/p0,  and-  p/p0  f°r  air. 
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Figure  3>_  Van   der  Waals '  equation  in  nondimensional  form. 
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Figure  5-_  Variation  of  the  ratio  of  specific  heats     y     for  a  Van  der 

Waals  gas.      TQ   =  520°  F  ats. 
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